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0\ : Abstract 
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f->i ; We show that neutrino oscillations in matter are always accompanied by collective 

' , forces on the medium. This effect may produce interesting consequences for the 

Q I background and the neutrino oscillations themselves. The force is maximal in the 

T^ • case of resonant neutrino conversion in the adiabatic regime. We study here the 

forces driven by I'e — '^fj,,T and Ue — Vg MSW conversion and shortly discuss their 

k><( ; possible relevance for the dynamics of a type II supernova. 
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It is well known that the neutrino propagation in matter can be drastically altered 
if neutrinos carry non-standard properties like non-vanishing mass, mixing, magnetic 
moment, and new interactions with the matter constituents and the classical fields. 
The MSW neutrino resonant conversion |^ and the resonant spin-fiavour precession 
(RSFP)P] are examples of phenomena that can take place due to the interplay of the 
coherent interaction of neutrinos with the matter and the anomalous neutrino properties. 
Usually, the study of such phenomena is focused on the fate of the neutrinos disregarding 
any possible implications for the medium itself. However, such an attitude could be not 
always justified especially in those cases where the neutrino fiux is very intense, like 
during a type II supernova (SN) explosion 0. 

Recent studies [^, Q showed that the coherent interaction of neutrinos with the mat- 
ter constituents can gives rise to macroscopic forces even in the case in which neutrinos 
are massless and they have only standard interactions. Such forces are usually called 
ponderomotive forces in analogy to the forces that arise in an electron gas in the presence 
of a non- uniform radiation field ||^, |^. A similar phenomenon takes place for the elec- 
trons, and other weakly interacting matter constituents, in the presence of nonuniform 
flux of neutrinos. According to the approach of Hardy and Melrose 0] , the background 
of neutrinos gives rise to a space-dependent self-energy, hence to a space-dependent mass 
correction, of the matter particles. By interpreting the statistical average of the mass 
correction as an interaction energy density U, Hardy and Melrose deduced the existence 
of a ponderomotive force per unit volume JF = — VW which is proportional to the Fermi 
constant Gp- The reader should not confuse the ponderomotive force with the more con- 
ventional force produced on the matter constituents by the incoherent elastic scattering 
off the neutrinos, which is proportional to G'jp. In the past years some work has been 
done concerning collective forces produced by coherent neutrino scattering in relation 
to possible detection of relic neutrinos f^. On the basis of the Born approximation and 
the assumption of spatially homogeneous neutrino flux, it was proved that to flrst order 
in Gp the effect is vanishing ^j. We note, however, that these arguments do not apply 
to objects of astrophysical sizes where the neutrino flux is not homogeneous. 

The relevance of neutrino induced ponderomotive forces for the physics of type II 
supernovae and other extreme astrophysical objects was investigated by several authors 
|§|, ^, H, |lOl- In some cases, non- negligible effects were predicted. 

The aim of this Letter is to show that neutrino oscillations in matter are always 
accompanied by collective forces on the background medium. Furthermore, we will show 
that in some phenomenological interesting cases, these forces are quite larger than the 
ponderomotive forces generated in the absence of neutrino oscillations. After deriving 
the general expression of the force produced by neutrino oscillations we shall focus on 
the case of MSW neutrino conversion. As an application, we shall briefly discuss the 
possible relevance of our results for the dynamics of a type II SN. 

Before entering into the details of our derivation, we have to mention that in the 
last few years some debate was going on concerning the correct expression of the pon- 



deromotive force in the absence of neutrino oscillations |TT| . Since our results might be 



only marginally affected by the outcome of this controversy, we do not enter here into 



such discussion. We observe only that our general treatment lead to results that are 
consistent with those of Hardy and Melrose. 

Similarly to what done by other authors, we determine the ponderomotive force per 
unit volume by taking the gradient of the interaction energy density U. The way we 
compute U is however more straightforward than that followed in Ref.[^ and it allows 
us to apply the formalism which is commonly used to treat neutrino oscillations in 
matter. For the sake of clarity we assume here that the relevant components of the 
medium are well described by some equilibrium distribution function. We also assume 
that at the position x, where the force is computed, neutrinos are already thermally 
decoupled from the background. According to the Liouville theorem neutrinos can still 
be described by a Fermi-Dirac distribution function fy^ (x; k) (where the variable x 
enters as a parameter) with an effective temperature which will depend on the distance 
from the decoupling position. The interaction energy density of the neutrinos with the 
i — th medium component is computed by taking the thermal average of the interaction 
Hamiltonian over the fermion distributions. We have 

— --/.„(x; k)P.„.,(x, k)K,,,(x, k) . (1) 



The quantity 



K„. = {^^\n,J^,)T (2) 

is the i/fl-potential in a medium composed only by the i-th matter component. In the 
above Ti^^i is the interaction Hamiltonian density and (\I/j|...|\I/j)t stands for the thermal 
average over the i-th particle component of the plasma. It should be noted by the reader 
that, in order to account for neutrino oscillations, in the equation (|I]) the integral over 
the neutrino momentum contains the probability P„^y^{'x,'k) = P(t'a —>■ t'b)(x, k). As 
we assume isotropic distributions for the neutrinos and the matter components and we 
disregard possible effects due to the polarization of the medium, it is easy to verify that 
in the one-loop approximation Vi^^i does not depend on the neutrino momentum k. In 
this approximation we find that the force per unit volume produced by all neutrino 
species on the i-th matter component is 



J-,(x) = -VW,(x) 

= JT(x)- E /7^P^a..(x,k)V[f.,(x;k)V.,i(x)] 



8L,b=e,fj.,T 

where 



■d^k 



•^r(x) = - E K,.(x)y— -^/.„(x;k)VP.„.,(x,k). (3) 

a,b=e,^,T V / 

The latter contribution JF°^"^ is the new term induced by the oscillations, or conversions, 
of a neutrino species into another - active or sterile. Clearly, that expression for J-'°^'^ in 
Eq. (^ accounts for whatever kind of neutrino conversion once the proper probability 
is specified. We note that in the absence of neutrino oscillations our expression (|) 



reproduces the results of Hardy and Melrose H, 11 



We predict several interesting effects of tlie ponderomotive force induced by neutrino 
oscillations. First of all, this force will produce a rearrangement in the density of the 
different matter components. Such an effect will be maximal in the case of MSW (or 
RSFP) neutrino resonant conversion. Indeed, in the resonance layer the neutrino survival 
probability undergoes a rapid variation which may give rise to strong macroscopic forces. 
For this reason, this is the case that we are going to discuss in more details here. As a 
secondary effect, the density modulation produced by the force will affect the neutrino 
propagation and the neutrino oscillations/conversion. We shall not study this second 
order effect here. However, it is worthwhile to observe that as a consequence of this 
effect and of its feedback on the ponderomotive force, non-linear effects appear which 



may lead to a sizeable energy transfer from the neutrino to the plasma |]T0| and to the 
amplification (or the fast damping) of the neutrino oscillations. 

1. Ponderomotive force due to MSW resonant conversion. It is worth recalling the 
main features of the resonant neutrino conversion. We consider the system of two neu- 
trinos z/g and Ux {x = fi,T,s) (here Ug denote an iso-singlet state), characterized by the 
difference of the eigenstate mass squares Sm^ = ml — ml and the vacuum mixing angle 
6. The coherent neutrino scattering off matter constituents can be described in terms 
of matter potentials as shown in Eq. (^. In the rest frame of the unpolarized matter, 
they read as: 

K.eT = ±^^(4sin2^^ + l)pY,^, (4) 

K,.eT = ±^^(4sin2^^ - l)py;^, (5) 

K.p = ^^(1 - 4 sin^ e^)pYp , (6) 

V2m„ 

Kan = f/ pYn, K,i = 0, i = e^,p,n (7) 

V2m„ 

here p is the matter density, Yi the concentration of the i-type component and m„ the 
nucleon mass (for the anti- neutrinos, the above potentials change sign). The Hamiltonian 
H governing the evolution equations in matter can be written as: 

_ ( ^-i^cos29-AV ^sin2^\ 
I ^sin2^ ; ' 



Grp 



. m„ 



'eff 



r) 



here A^ is the effective matter potential for the system of z/g — Vx and pes = pY where 
the net effective concentration F is F = Fg = Fg- — Yg+ for the fg — v^^^ system and 
Y = Ye — \Yn for the p^ — ^s channel. The mixing angle 6m and the neutrino wavelength 
Am in matter are given by: 

3i^2 2^ ^ {5m-^u.2ef 

"^ (5m2 cos 29 - 2EAVY + {5m^ sin 29f ' ^ ' 



A„. = ^"'^ (10) 

'{5m'^ cos 26 - 2E^Vf + {5m'^ sin 2^)2 



where A = -j — - is the vacuum wavelength. 

As it is well known the efficiency of the conversion is determined by the resonance 

. ^m^ cos 26* . ^^ - - , ,. , . . d6'in ,, , 

condition, = AV and by the adiabaticity property — — <^ vr/Am, where 

2E dr 

^ = sm=2«)„^4^ft-'(r). (11) 

dr om^ sm 26 

;,-!(,) ^ d Injpe.) ^ 

dr 
The t^e survival probability at a certain distance r from the source is given by 

P,Xr\ ^) = ^ + (^ - Pc) cos26Uri) cos2^„(r) , (13) 

where Tj is the (radial) coordinate of the neutrino source and the function P^ is the 



probability of jumping from one matter eigenstate to the other |]13|. In typical cases, 
EAV ^ 6m'^ in the neutrino production region which implies cos26m{ri) ~ —1. We 
assume the resonant conversion to be adiabatic -Pc — (this will be justified later on). 
Then, from Eqs. (|rT]) and (|T^) we can easily obtain 

^ = sin2^^^ = sin^ 26^ , T^ ^, h-\r) . (14) 

or dr 0171'^ sm 26 



We note that the r.h.s. of Eq.(|l4D attains a maximum when the resonance condition is 
fuimied (sin2 2^^ = l). 

Now we are ready to turn back to the expression for the ponderomotive force induced 
by the neutrino resonant conversion given in Eq. (^. For simplicity, we assume that 
only z/g's are created by some point-like source with a thermal distribution fu^lr; E) and 
propagate isotropically trough an inhomogeneous medium like that e.g. present in a 

star. Expressing the neutrino density in term of the luminosity Li, {ny = -^^ — -) we 

ATcr^{E) 

find from Eq. (^ 

L IdEE'U{E)^^^^p^ 
A W [v,^, '^^-'^47rr2 JdEE^U^{E) 

= ~{Vy.^-V,J-^^^^^h-\r)I{r) (15) 

dm? sm 26 Airr^ 

where 

■<'') " JlEE^KW) ■ ''«> 



I/tan2fl 




10"-' 10" 



Figure 1: The ratio /*(6')/tan2^ is plotted as a function of sin^ 26. 



Since the function I is obtained by a convolution of siv? 26m{E) with the neutrino dis- 
tribution function, it attains the maximum value I* at the position r^, = r.^es\E={E)i 
where {E) is the neutrino mean energy. The actual value of I* depends on the width 

dTTL 

AE oc sin 29 of the function sin^ 29m{E), hence on the degree of adiabaticity of the 
neutrino conversion. Therefore the maximal force is 



^rw = -(K..-K..)- 



Inr'^iE) 



h 



-1/ 



, h{9) 
tan26' 



(17) 



By comparing this expression with that of the conventional force driven by the gradient 
of the neutrino density (i.e. the second term in eq. (^, hereafter named JF^^^) we note 
a significative difference. While the latter is proportional to the net neutrino density A^,^, 
the former depends only on the helicity species which undergoes the resonant conversion. 
Such a different behavior may play a crucial role in SNs and in other astrophysical 
or cosmological environments. Upon making the comparison more quantitatively, we 
consider the forces produced on the electron component of the medium. We find 



J^° 



jrHM 

{E) 



nh-^ h.{9) 
4^ tan2e 



where the parameter ^ = 1 — ^ accounts for the possible difference in the u^. and z/g 
spectra (we assume that Ly^ = LpJ. In the Fig.l we plot the function /^,(^)/tan26'. 
We observe that this function does not lead to any suppression in the phenomenological 
interesting range 10~^ < sin^ 29 < 2 x 10^^. The actual value of r^h~~^ depends on the 
matter density profile and on the kind of neutrino transition (active-active or active- 
sterile) considered. Below we will show as in a SN such a quantity typically ranges from 
a values of few units to some powers of 10. We should keep in mind, however, that 
the adiabaticity requirement implies an upper limit to the allowed value of h'^. This is 
given by 



h-^ < h-^ 



adiab 



5m^ sin^ 29 
2{E)cos9 



(19) 



Above this value the scale setting the gradient of P^^ is Am which, in the non-adiabatic 

regime, is larger than h. Therefore the maximal force is achieved in the adiabatic regime. 

For completeness, we also compare the force J^"^'^ with the outward force due to the 

q2 ^ 

elastic scattering off electrons, JF^^^^** ~ — - — '^iE)Tene. We find 

Ff ^t\^ '^ \^ Te P2 l0-5cm-i; tan2^ ^ ' 

where Te is the electron temperature. Clearly J^°^^ ^ J^^^^^ . We observe, however that 
these forces are qualitatively quite different and cannot be always directly compared. 
For example, the ponderomotive force naturally excite acoustic and plasma waves which 
can hardly be done by the force due to the incoherent neutrino scattering. Therefore 
the former may induce effects which are not produced by the latter allowing its possible 
identification. 

2. Ponderomotive force in supernova. As an application of our previous results, we 
now investigate what kind of effects the ponderomotive forces produced by neutrino 
resonant conversion may give rise to during a type II SN explosion. We consider the 
region above the neutrino sphere at a time after the core bounce. In this epoch all 
the neutrino (as well anti-neutrinos) species are emitted with approximately the same 
luminosity. However the individual neutrino energy distributions may be quite different. 
For this reason the force produced, say by the u^ — > v^ MSW conversion is not canceled 
by the opposite force produced by v^ -^ v^. that will occur more externally. In both cases 
the resonance takes place in a dynamically relevant region for brr? = 10 ^ 10^ eV^ if 
(E) ^ 10 MeV. In the earlier epoch (t < 1 s) with a typical profile |T^ A^e = lO^^rf ^ cm"'^ 



[r-j = r/10^ cm), from eq.([l8|) we find 



rose 

e 



1 with ^ ~ 0.3. At later times, as the 



JTHM 

scale h becomes smaller, this ratio can increase up to an order of magnitude. 

A different picture may emerge for the z/g — z/g, or the u^ — i^s, channels. We recall that 
in this case pefr = p(3Ve ~ 1) and the resonance condition can be satisfied for arbitrarily 
small values of Sni^. Hence, for small values of 6m'^ the resonance take place nearby 
the point where the potential approaches zero. As a consequence /i~^(r*) could be quite 
large. At the same time, however, the adiabaticity requirement (plQ]) set an upper limit 
to h~^ which is proportional to Sm"^. Adopting standard SN density profile |jl5[ and 
assuming sin^ 29 < 10^^ , we found that the maximal force is obtained for Srn^ ~ 1 eV^, 
corresponding to AV ~ 10"'^ eV and h~^ ~ 10~^ cm~^ (r* ~ 150 km). Therefore the 
ratio (ITB|) can be as large as 10^. 



A comment at this point is in order. Because the electron neutrino potential changes 
sign at the radius where Ye = 1/3, two kinds of resonance should take place close to 
this position. If we assume < Srn^ < 10^ eV^, and since Y^ is a growing function of 
the radius above the neutrino-sphere, the transition u^ — i^s will first occur in the region 
where Ye < 1/3. It is easy to verify that in this case the force pulls electrons inwards. 
Soon afterwards the transition z/g — z/^ take place in the region with Ye > 1/3 and thereby 
the force goes in the opposite direction. It happens |T5[ that the two types of conversion 



take place very close to each other. However the corresponding forces should not cancel 
as the average energy of v^ and Ve. differs of some 30-50%. In the case bvr? > 10^ eV^ 
only the conversion z^e — Vs can happen as far as Ye < 1/3 |T5[| . 



It is interesting to estimate the velocity which background matter close to the res- 
onance position acquires under the action of T^^ in a time interval of the order of the 
dynamical time ~ 0.1 s. A straightforward computations shows that velocities as large as 
~ 10^ km/s can be reached. Other interesting effects arise if these ponderomotive forces 
induce plasma instabilities |jlO| which may give rise to a significative energy transfer 



from the neutrino to the plasma. 

In conclusion, we think that the phenomenological perspectives open by the study of 
the ponderomotive force produced by neutrino oscillations are rich and deserve further 
investigation. 
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